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ABSTRACT 

We have used the Berkeley-Illinois-Maryland-Association (BIMA) millimeter array 
outfitted with sensitive cm-wave receivers to search for Cosmic Microwave Background 
(CMB) anisotropics on arcminute scales. The interferometer was placed in a compact 
configuration which produces high brightness sensitivity, while providing discrimination 
against point sources. Operating at a frequency of 28.5 GHz, the fwhm primary beam 
of the instrument is ~ 6.6'. We have made sensitive images of seven fields, five of which 
where chosen specifically to have low IR dust contrast and be free of bright radio sources. 
Additional observations with the Owens Valley Radio Observatory (OVRO) millimeter 
array were used to assist in the location and removal of radio point sources. Applying 
a Bayesian analysis to the raw visibility data, we place limits on CMB anisotropy flat- 
band power Q flat = 5.6±|;g /iK and Q f i at < 14.1 /xK at 68% and 95% confidence. The 
sensitivity of this experiment to flat band power peaks at a multipole of £ = 5470 which 
corresponds to an angular scale of ~ 2 . The most likely value of Qfiat is similar to the 
level of the expected secondary anisotropies. 
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1. Introduction 

The Cosmic Microwave Background (CMB) has 
the potential to be a powerful probe of the early uni- 
verse. In the standard inflationary model, the CMB 
is imprinted with anisotropics that reflect the dis- 
tribution of matter at the epoch of recombination. 
Observations of the CMB at degree angular scales 
probe structures which have recently collapsed at this 
epoch and for which the distribution of anisotropics 
is extremely sensitive to the cosmological model. On 
scales smaller than a few arcminutes, photon diffu- 
sion and the finite time of recombination damp the 
primordial fluctuations to near zero amplitude (Hu 
Ic White 1997). However, the subsequent reioniza- 
tion of the universe can create a host of secondary 
anisotropies of the CMB. For a review of the subject 
see Haiman & Knox (1999). On arcminute scales, 
secondary anisotropies generated since recombination 
are likely to dominate the primary signal. 

There have been many previous searches for an- 
isotropy in the CMB on arcminute scales; both inter- 
ferometric and single dish techniques have been used 
successfully. In this paper we describe a search for 
arcminute-scale CMB anisotropies with the Berkeley 
Illinois Maryland Association (BIMA) interferometer 
in a compact configuration at 28.5 GHz. We begin 
with a discussion of the instrument, field selection 
and observations in §||. In §[| initial data reduction is 
described including point source detection and mea- 
surement. The Bayesian maximum likelihood analysis 
we apply to the data is described in §||. The results of 
applying this formalism to our data are presented in 
§H including a discussion of the effects of point source 
subtraction. In §j(| we discuss the levels of the ex- 
pected signals. We summarize the previous work in 
this field in Finally, in we summarize the re- 
sults and discuss prospects for future observations. 

2. Instrument & Observations 

The advent of low-noise, broad-band, millimeter- 
wave amplifiers has made interferometry a particu- 
larly attractive technique for detecting and imaging 
low contrast emission, such as anisotropy in the CMB. 
An interferometer directly samples the Fourier trans- 
form of the intensity distribution on the sky. By 
transforming the interferometer output, images of the 
sky are obtained which include angular scales deter- 
mined by the size and spacing of the individual array 
elements. In this section, we describe the BIMA in- 



strument, the selection of the fields, and their obser- 
vation. 

2.1. Instrument 

The anisotropy observations described here were 
made with nine elements of the BIMA array out- 
fitted with sensitive cm-wavelength receivers. The 
BIMA antennae are 6.1m in diameter and produce 
6.6' beams at 28.5 GHz. The receivers are based on 



low noise HEMT amplifiers ( Pospieszalski et al. 1998) 
and arc configured to respond only to right circularly 
polarized radiation. At the operating frequency of 
28.5 GHz, the receiver noise temperatures range from 
13 — 18 K and the system temperatures were typi- 
cally 35 — 55 K depending on source elevation. The 
signals from the individual array elements are com- 
bined in the BIMA 2-bit digital correlator which was 
configured for these observations with 8 contiguous 
100 MHz sections of 32 channels each. The array el- 
ements were placed in a compact 2-dimensional con- 
figuration which provides high brightness sensitivity 
as well as sufficient resolution to identify radio point 
sources. 

We have also made supporting observations with 
the Owens Valley Radio Observatory (OVRO) array 
outfitted with the same cm-wave receivers. These 
observations were used primarily for locating and 
measuring point sources in the observed fields. The 
OVRO antennae are 10.4 m in diameter and produce 
3.8' beams at 28.5 GHz. For a description of the 
instrument see Carlstrom et al. (1996). The BIMA 
and OVRO arrays, outfitted with cm- wavelength re- 
ceivers, have been used to image more than 20 clusters 
(|Carlstrom et al. 1999 



2.2. Field Selection 

In this paper, we report the results of a search for 
CMB anisotropy in seven independent fields. Two 
of the fields, PC1643+46 and VLA1312+32, are cen- 
tered at previously reported microwave decrements 



(lones et al. 1997; Richards et al. 1997). In a compan- 
ion paper (Holzapfcl ct al. 1999), we demonstrate that 
our data are inconsistent with models used to describe 
the claimed decrements. In 1997, we observed a third 
field in the direction of the quasar PSS0030+17 which 
was originally selected as a distant cluster candidate. 
The quasar is at redshift z — 4.32 and ha s two Lya 
break galaxies within 10" ( Djorgovski 1997| ). The way 
these three fields were selected would prevent us from 
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claiming that they could produce an unbiased mea- 
surement of CMB anisotropy. However, for the pur- 
pose of placing upper limits on CMB anisotropy, we 
are justified in making use of these observations. 

In 1998, four additional fields were selected to be 
evenly distributed in right ascension and at conve- 
nient declinations for observations with the BIMA ar- 
ray. The fields were chosen to be in regions of low dust 
emission and contrast as determined from examina- 



of Mars which is uncertain by approximately 4% (see 
discussion in Grego 1998). 



tio n of IRA S 100 ^m map s. The VLA NVSS QCondon 
et frl. 1998| ) and FIRST flWhite et al. 1997| ) surveys 
were then used to select regions free of bright point 
sources at 1.4 GHz. Finally, we used the SkyView^J 
Digitized sky survey and ROSAT WFC maps to check 
for bright optical or x-ray emission which could com- 
plicate follow-up observations. The pointing centers 
for each of the 7 fields are given in Table [|. 

For one of the new fields, we used the OVRO ar- 
ray to check for radio sources at 28.5 GHz before be- 
ginning observations with the BIMA array. With 
a 7 pointing mosaic, we reached a map rms flux 
of ~ 120 //Jy over a 8' region containing the entire 
BF0028+28 field observed with BIMA. We discovered 
a single source with a flux of 1.4 mJy. This and all 
fluxes in this paper have been corrected for the atten- 
uation of the primary beam unless otherwise specified. 
The pointing center for the BIMA observations was 
chosen so that this source lied outside the observed 
field of view. Unfortunately, we did not have time to 
image all the observed fields with the OVRO array 
and it is possible that some of the fields observed in 
1998 suffer from low level point source contamination. 

2.3. BIMA Observations 

All observations were made during the summers of 
1997 and 1998, interspersed between observations of 
the Sunyaev-Zel'dovich effect (SZE) in x-ray selected 
clusters. In 1998, we selected four fields spaced in 
right ascension so that at any given time one of them 
had a hour angle suitable for observation. Each 20 
minute source observation was bracketed by the ob- 
servation of a calibration source. Including the time 
for calibration cycles, the fraction of time spent on 
source was ~ 60%. The integration times for each 
of the 7 fields are given in Table § The fluxes of 
the calibration sources are all referenced to the flux 
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(pttp:/ /skyvicw.gsfcnasa.gcrv) located at NASA Goddard 



Space Flight Center. 



3. Editing and Inspection 

The data are edited using several criteria to ensure 
the integrity of the calibration and that the results re- 
main free of systematics. It is possible for the beam of 
one dish to be obscured by one of its neighbors in the 
array. Baselines involving telescopes within 3% of the 
shadowing limit are discarded. The spectral channels 
are inspected for any interference and removed if they 
are believed to be contaminated. Low signal to noise 
channels near the edges of the correlated bandwidth 
are not used. The effective noise bandwidth of the 
correlator after accounting for the 2-bit digitization 
and removed end channels is ~ 540 MHz. Records 
with spurious system temperatures, caused by failed 
or aborted calibration cycles are discarded. Source 
data not bracketed by successful calibration cycles are 
discarded. During periods of poor weather, the phase 
coherence of the calibration sources becomes poor. 
All data that are bracketed by calibration cycles with 
poor phase coherence are also discarded. 

3.1. Point Sources 

Each data set is transformed to create a map with 
the DIFMAP package (Shepard, Pearson, & Taylor 
1994). The maps are then searched for statistically 
significant unresolved emission. In order to remain 
unbiased in our search for point sources, we limit our- 
selves to the range of baselines greater than 2.4 kA 
(£ > 15000), which are completely independent of the 
baselines used in the anisotropy analysis. That way, 
we can be assured that the anisotropy results will not 
systematically depend on the point source detection 
and subtraction. In general, we find the flux and po- 
sitions of the sources by fitting the Fourier transform 
of source model directly to the visibility data. The 
source model is multiplied by the measured primary 
beam response to take the attenuation of the source 
into account. In Table ||, we list the positions and 
fluxes of the significant point sources in the observed 
fields. 

The brightest point source discovered was centered 
at a = 00 h 30 m 37.0 s , S = +17° 05' 12" (J2000), offset 
from the pointing center of the PSS0030+17 field by 
Aa = +294" and A5 = +152". The observed flux 
is attenuated by the finite size of the array element 
beams which are measured to have a fwhm of ~ 396 . 
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Field Positions and Observation Times 



Fields a (J2000) 5 (J2000) Observation Year(s) Time (Hours) 



PC1643+46 16 h 45 m 11.3 s +46° 24' 56" 1997, 1998 43T 

VLA1312+32 13 h 12 m 17.4 s +42° 38' 05" 1997 35.5 

PSS0030+17 00 h 30 m 16.4 s +17° 02' 40" 1997 36.6 

BF0028+28 00 h 28 m 04.4 s +28° 23' 06" 1998 77.6 

HDF1236+62 12 h 36 m 49.4 s +62° 12' 58" 1998 32.0 

BF1821+59 18 h 21 m 00.0 s +59° 15' 00" 1998 43.5 

BF0658+55 06 h 58 m 45.0 s +55° 17' 00" 1998 44.7 



Table 1: Coordinates of the observed fields, years of observation, and cumulative time on source. 



Subtracted Point Sources 







Positions (J2000) 


From Map Center 




Field 


Source 


a 


5 


Aa 


AS 


Flux (/zJy) 


PSS0030+17 




00 h 30 m 37.0 s 


+17° 05' 12" 


294 


152 


12800 ± 1600 


HDF1236+62 




12 h 36 m 44.4 s 


+62° 11' 33" 


-34 


-85 


347 ± 59 


PC1643+46 


S[ 


16 h 45 m 21.1 s 


+46° 25' 46" 


102 


50 


227 ± 70 


PC1643+46 


S! 2 


16 h 45 m 12.2 s 


+46° 24' 08" 


9 


-48 


113 ±62 


PC1643+46 


S' 3 


16 h 45 m 22.7 s 


+46° 24' 17" 


118 


-40 


237 ± 72 


PC1643+46 


Si 


16 h 45 m 20.9 s 


+46° 25' 43" 


100 


47 


227 ±63 


PC1643+46 


s 2 


16 h 45 m 12.5 s 


+46° 24' 13" 


6 


-43 


222 ± 63 


PC1643+46 


s 3 


16 h 45 m 23.0 s 


+46° 24' 21" 


115 


-36 


262 ±71 



Table 2: Coordinates, distances from the anisotropy map center, and intrinsic fluxes of the significant sources in 
the field at 28.5 GHz. For PC1643±46, we first show the results when the positions of the sources (primed) are 
fixed and second when the positions are allowed to assume their best fit values. 



This source is far from the center of the map and 
has an observed flux of 1.5 mJy. Correcting for the 
primary beam response, the intrinsic source flux is 
determined to be 12.8 ± 1.6 mJy. 

A significant point source was detected in deep ob- 
servations of the HDF1236±62 field with both the 
OVRO and BIMA arrays. There were two pointings 
of the OVRO array with one centered on the posi- 
tion of the suspected point source. We simultaneously 
fit the raw OVRO and BIMA visibility data with a 
single component point source model. The source 
is the brightest in the observed field with a flux of 
345 ± 59 fiJy. It was centered at a = 12 h 36 m 44.4 s , 
S = +62° 11' 33" (J2000), offset from the BIMA 
pointing center by Aa = —34 and AS = —85 . This 
source has the same position as the brightest source 
found in a deep radio image of the HDF at 8.4 GHz 



The field PC1643±46 was pre viously imaged wit h 
the Ryle telescope at 15 GHz ( Jones et al. 1997 ). 



(Richards et al. 1998) 



Three point sources were found in the observed field. 
In addition to our observations with the BIMA ar- 
ray, we have also imaged this same region with three 
pointings of the OVRO array operating at 28.5 GHz. 
Each of the three pointings were chosen to place the 
map center near one of the suspected point sources. 
The field observed with BIMA was imaged with rms 
flux density ranging from 40 — 100 /xJy beam -1 . We 
have performed a simultaneous fit to the BIMA and 
OVRO visibility data in order to determine the posi- 
tions and fluxes of the sources at 28.5 GHz. We have 
determined the fluxes of the sources using two dif- 
ferent methods. First we fixed the positions of the 
sources to the positions found with the Ryle obser- 
vations and solved for the three source fluxes. In- 
troducing these three free parameters decreased the 
X 2 of the fit to the visibility data by 25; this indi- 
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cates that the field suffers significant contamination 
from these sources. The uncertainties for the source 
fluxes, shown in Table correspond to the change in 
flux which results in a increase in \ 2 of one while the 
other free parameters (two other fluxes) are allowed 
to assume their best fit values. 

We repeated this analysis allowing the source po- 
sitions as well as fluxes to vary. This procedure al- 
lows for differences in the positions determined by the 
BIMA and OVRO analysis and those of Jones et al. 
(1997). By allowing the positions to vary, and adding 
6 new free parameters to the model, the \ 2 of the fit 
is reduced by 14 from the case where the source po- 
sitions are fixed to the Ryle positions. Therefore, the 
differences from the Ryle positions are significant and 
in the rest of this work we adopt these new source po- 
sitions. The differences between each of the new po- 
sitions and those of Jones et al. (1997) are less than 
6". When the uncertainties in the Ryle positions and 
those found here are taken into account, the positions 
determined by the two experiments are found to be 
consistent. The errors for each of the point source 
fluxes correspond to the change in flux required to 
produce to a change in \ 2 of one while the free para- 
meters (flux of the other two sources and the posi- 
tions of all three sources) are allowed to assume their 
best fit values. In Table ||, we list the measured po- 
sitions and fluxes of all the detected sources. These 
sources are removed from the raw data by taking the 
Fourier transform of the point source model modu- 
lated by the primary beam response and subtracting 
it directly from the visibility data. 

3.2. Image Statistics 

We have produced and analyzed images for each 
of the observed fields. The results for the long base- 
line data used in point source subtraction are listed 
in Table ||. We limit the data to baselines > 2.4 kA to 
guarantee that the data used to determine the point 
source fluxes and positions are completely indepen- 
dent of the anisotropy data. The rms for all the data 
is considerably lower. The map rms indicates the ac- 
curacy with which the flux of point sources can be 
measured with this subset of the BIMA data. For the 
fields PC1643+46 and HDF1236+62 the point source 
sensitivity is considerably better than listed here due 
to the supporting OVRO observations. The results 
using only the short baselines used in the anisotropy 
analysis are listed in Table |]. For the short baseline 
maps, we also express our results in terms of the rms 



Rayleigh- Jeans (RJ) temperature fluctuations. For 
both the short and long baseline maps, the observed 
rms values are compared to those expected from the 
noise properties of the visibilities. For the short base- 
line results, there are approximately ten independent 
beams in each observed field. If we assume that the 
map values are dominated by Gaussian distributed 
noise, the measured rms should be the same as the 
estimated value within approximately 10% at 68% 
confidence. 

4. Analysis 

Several recent papers have dealt with the anal- 



ysis of CMB data from interferometers (Martin & 


Partridge 1988 


; Subrahmanyan et al. 199? ; Hobson, 


Lascnby & Jones 1995; Hobson & Mageuijo 1996; 


Par- 


tridge et al. 1997; 


White et al. 199£). In this work, we 



follow the formalism presented in White et al. (1998) 
for the Bayesian analysis of CMB data. 

In theories which predict Gaussian temperature 
fluctuations, the fundamental theoretical construct is 
the correlation matrix of the measured data. Since 
the data are the visibilities measured at a set of points 
Uj, we will need to know the correlation matrices for 
the signal and noise of the observed visibilities. The 
measured fluxes are given by 

r) Ft f 

-Q^Jdx AT(x) A(x)e 2 — x , (1) 



V(u) 



where AT(x) is the temperature distribution on the 
sky, A(x) is the primary beam of the telescope, 



9B„ OI fk B T\ 2 x A e x 
2k B 



dT 



he 



(e» - l) 2 



(2) 



k B is Boltzmann's constant, and x = hv /fcsT cm b. We 
define the visibility correlation matrix, 



C% = (V*( Ui )V( Uj )) 



(dB v 



\ dT 



wdw S{w)Wij(w), (3) 



which is proportional to the product of the power 
spectrum, S(w), and the visibility window function. 
The window function is given by 

p2rr 

Wij(\w\)= d0 w A*(u i -w)A(u j -w) , (4) 
Jo 

where A(u) is the Fourier transform of the telescope 
primary beam. In the case of a single flat band power 



5 



Image Analysis for Baselines > 2.4 kA 



Field 


beamsize( ) 


rms (fiJyh 


earn ) 










estimated 


measured 


PC 1643+46 


18.9 


X 


25.0 


132 


139 


VLA 1312+32 


18.8 


X 


26.0 


134 


122 


PSS 0030+17 


18.1 


X 


25.2 


155 


150 


BF 0028+28 


21.2 


X 


21.6 


136 


127 


HDF 1236+62 


22.5 


X 


23.9 


195 


207 


BF 1821+59 


21.7 


X 


23.6 


152 


160 


BF 0658+55 


21.6 


X 


23.2 


158 


157 



Table 3: Image statistics for maps created using only the long baselines used to measure point sources. Columns 
3 & 4 list the estimated and observed map rms. 



Image Analysis for u-v Range 0.63 — 1.2 kA 

Field beamsize( ) rms (/iJy beam -1 ) rms (/iK) 

estimated measured estimated measured 



PC 1643+46 


98.3 x 116.1 


187 


191 


24.6 


25.1 


VLA 1312+32 


95.2 x 113.4 


225 


223 


31.3 


31.0 


PSS 0030+17 


99.1 x 115.5 


202 


203 


26.5 


26.6 


BF 0028+28 


108.9 x 118.8 


127 


99 


14.7 


11.5 


HDF 1236+62 


110.9 x 122.0 


206 


168 


22.8 


18.6 


BF 1821+59 


108.4 x 122.4 


174 


224 


19.7 


22.4 


BF 0658+55 


105.2 x 148.0 


252 


304 


24.3 


29.3 



Table 4: Image statistics for maps created using only the short baselines used in the anisotropy analysis. Columns 
5 & 6 give the image results in rms RJ map temperature. 



and £ > 60, we can write 
fOD v \ 



5tt \0T J 



rw 



WjAw\ ■ [5] 



where 



2/^^V(M I) 



(6) 



is the normalization of the power spectrum. The cor- 
relation function of the noise is diagonal with elements 
given by 



1 



(7) 



where <7i is the variance of the measured visibilities. 
For a given set of n measured visibilities, one can 
test any theory, or set of {Ce}, by constructing the 
likelihood function (for complex visibilities) 

c{{Cf}) = ^htc exp [-VM^vM] , 

(8) 



where C lj = CY j 



C£f is the correlation matrix of 



19951 ) 



sibilities at and (Hobson, Lasenby fc Jones 



4.1. Joint Confidence Intervals 

Invoking Bayes' theorem and assuming a uniform 
prior for the amplitude of the fluctuations, we can 
determine the probability that the correct result is 
contained in an interval /, 



P{I) = 



JjC(z)dz 
J °° C(z)dz 



(9) 



The confidence interval corresponding to a probability 
Pq is given by the Iq such that P(Iq) = Po and C[z € 
lo] > £[ z ^ ^o]- If the fields are entirely independent, 
the joint likelihood for the combination of the data 
sets to be described by a given model is simply equal 
to the product of the likelihoods for the individual 
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data sets, 



C(z) = HCi(z) . 



(10) 



Combining the diagonal window functions correspond- 
ing to each visibility weighted by the noise, we can 
construct an effective diagonal window function to de- 
termine where the experiment is most sensitive; 



v W u {i) 
2^1 £ 2^ 



where Wi = l/of . With this normalization, 

= 



(ii) 



(12) 



Using the data weighted window function, we can de- 
termine the effective multipole of the experiment as- 
suming the power spectrum is flat; 



(13) 



4.2. Binning 



The data sets for each of our fields contain on the 
order of N = 10 5 visibilities. The inversion of the 
covariance matrix required for the analysis is a iV 2,8 



process (Press et al. 1996). As discussed in Hobson, 
Lasenby & Jones (1995), considerable compression of 
the data is necessary if the analysis is to be completed 
in a reasonable amount of computing time. We divide 
the u-v plane into a grid of cells; all the visibilities 
within a cell are combined, weighted by the reciprocal 
of their estimated noise variance; 



V a = 



E 



Vi 



E ■ 



(14) 



Wc determine noise weighted u-v positions just as we 
have determined the values for the visibilities; 



EUi \^ Vi 



E ^ E 



(15) 



The sampling theorem tells us that the Fourier trans- 
form of the sky is completely specified by a sampling 
of the u-v plane on a regular grid with Am, At; < 
l/29 p where 6 P is the angular radius of the primary 



beam. Because our final u-v points do not form a reg- 
ular grid, it is possible that we would have to decrease 
the size of a our grid by a factor of two to strictly 
satisfy this criterion. Following Hobson, Lasenby & 
Jones (1995), we define the extent of the beam as the 
point at which the beam falls to 1% of its peak value, 
dp = 7.45'. We determine that the u-v plane must be 
sampled more densely than Am, Aw = 230. In prac- 
tice, this is simple to achieve and we sample the u-v 
plane at intervals of Au, Av — 60 in the analysis pre- 
sented in this paper. This compresses the number of 
u-v points to ~ 200 for each of the data sets. 

4.3. Calibration 

The likelihood analysis code has been checked by 
the analysis of simulated data sets. We took a data 
file from one of our observations, removed the visi- 
bilities, and replaced them with Gaussian distributed 
noise with the same weights as the original data. To 
this we added the Fourier transform of a realization 
of CMB anisotropy with flat band power that had 
been modulated with the measured BIMA primary 
beam. We produced 100 such data sets each with 
independent CMB and noise realizations. The sim- 
ulated data were then analyzed in exactly the same 
way as the real data treating each data set as an in- 
dependent (uncorrelated) observation. For an input 
flat-band power with Qf\ a t — 30 /xK, we found that 
the likelihood peaked at Qfiat = 31 /xK with values 
29 - 33 fjK and 27 - 35 fxK allowed at 68% and 95% 
confidence. We interpret this as a demonstration that 
the analysis code is correctly calibrated. 

5. Results 

In this section, we use the method that we have 
described above to determine the relative likelihoods 
that the observed fields are described by a model for 
the CMB fluctuations with flat band power Qfiat- 

Before proceeding, we discuss the effect of sub- 
tracting the known point sources on the anisotropy 
results. Three of the observed fields are known to have 



significant point source contamination. In Section 3.1 



we discussed the determination of the source fluxes 
and positions. For the three contaminated fields, we 
subtract the Fourier transform of the point source 
model from the raw visibility data. In Figure [l], we 
plot the results of the likelihood analysis for the field 
PC1643+46 before and after the subtraction of the 
three detected point sources. The results are normal- 



7 



ized to unity likelihood for the case of no anisotropy 
signal. We have performed the same analysis with the 
point source fluxes found when the positions of the 
sources are fixed to the results of Jones et al. (1997). 
The results are essentially identical; the most likely 
value for Qfiat is again zero. In all further analysis, 
we use the point source fluxes and positions found 
from the fits to the OVRO and BIMA data, In Ta- 
ble [|, we show the results for Q fi at for the analysis 
of the three contaminated fields before point source 
subtraction. It is clear that neglecting to subtract the 
known point sources can lead to an erroneous detec- 
tion of anisotropy. 




Fig. 1. — The relative likelihood that the observed 
signal in the PC1643+46 field is described by flat 
band power with amplitude Qfiat- The dotted line 
corresponds to an analysis ignoring the measured 
point sources and the the solid line is the result when 
the measured point source fluxes are subtracted from 
the visibility data. 

In Figures || and ||, we plot the relative likelihood 
for each of the observed fields, where we have sub- 
tracted the detected point sources listed in Table 0. 
In Table [| we list the 68% and 95% confidence inter- 
vals in Qfiat for each of the observed fields. The fields 
are independent and we can apply eqn. |l^ to deter- 
mine the joint likelihood for the combination of fields. 
The relative likelihoods of the joint fits are plotted in 
b lgure § In Table |, we also list the confidence inter- 
vals in Qfiat found from the joint fits. Because of the 
different array configurations and declinations of the 
sources, the window function for each observation is 
slightly different. 



data. These window functions are plotted as a func- 
tion of multipole in Figure g. Finally, we have used 
eqn. 

4// 



13 to determine the effective multipole number, 
5470, of the results quoted here. 
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to determine 



We have used eqn. 
effective diagonal window functions corresponding to 
the 1997, 1998, and combination of the 1997 and 1998 



Fig. 2. — The relative likelihood that the observed 
signal in each of the 1997 data sets is described by flat 
band power with amplitude Qfiat- The solid, dashed, 
and dotted lines are for the fields PC 1643+46, 
PSS 0030+17, and VLA 1312+32 respectively. 

One is immediately struck by the fact that the joint 
likelihood for all the data, shown in Figure |J, peaks 
at Qfiat > 0. However, this result has fairly low sig- 
nificance. Each of the individual fields are consistent 
with no signal on the sky at 95% confidence and the 
confidence of a non-zero Qfiat for the joint likelihood 
is only 44%. 

After point source subtraction, none of the fields 
observed in 1997 produce a non-zero anisotropy sig- 
nal. This is further evidence that the PC1643+46 and 
VLA1312+32 data are inconsistent with the presense 
of massive galaxy clusters in these fields as has been 
demonstrated by Holzapfel ct al. (1999). There is also 
no evidence for a distant cluster in the PSS0030+17 
field. From Table g, we can see that including the 
1997 data in the joint likelihood only reduces the lim- 
its on Qfiat- 

5.1. Limits on Point Source Contamination 

All of the excess power is found in the two fields 
which have the poorest limits on point source contam- 
ination of any of the seven fields. It is possible that 
point source contamination contributes to the signal 
in these fields. If the observed power is dominated by 
Poisson distributed point sources, we expect to find 
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Results Before Point Source Subtraction 







Qflat (juK) 




Field 


Most Likely 


68% 


95% 


HDF 1236+62 


4.8 


0.0-17.8 


0.0-33.2 


PSS 0030+17 


11.8 


0.0-31.6 


0.0-61.8 


PC 1643+46 


18.2 


4.8-29.0 


0.0-42.2 



Table 5: Qfiat results before subtraction of known point sources. Compare with Table |6J to see the results after 
point source removal. 



Most Likely Q f\ a t and Confidence Intervals 







Qflat(lM) 






Field 


Most likely 


68% 


95% 


BF 0028+28 


0.0 


0.0 


- 8.5 


0.0- 


17.5 


HDF 1236+62 


0.0 


0.0 


- 13.0 


0.0- 


26.6 


BF 1821+59 


18.6 


8.8 


- 28.8 


0.0- 


37.8 


BF 0658+55 


19.0 


8.4 


- 30.4 


0.0- 


40.6 


Combined 1998 fields 


8.8 


2.4 


- 12.8 


0.0- 


17.4 


PC 1643+46 


0.0 


0.0 


- 13.1 


0.0- 


26.6 


VLA 1312+32 


0.0 


0.0 


- 20.5 


0.0- 


38.0 


PSS 0030+17 


0.0 


0.0 


- 14.4 


0.0- 


29.2 


Combined 1997 fields 


0.0 


0.0 


- 8.9 


0.0- 


17.4 


All fields 


5.6 


0.0 


- 9.6 


0.0- 


14.1 



Table 6: Results of the Bayesian analysis for each of the observed fields. Column 2 gives the most likely value for 
band power amplitude Qfiat- Columns 3 & 4 give the 68% and 95% confidence intervals. 



Q fiat <x We have verified this scaling by the anal- 
ysis of both real and simulated data with significant 
point source emission. If this scaling was observed in 
these data sets, it would be a "smoking gun" for point 
source contamination. 

Two of the fields we have observed, BF0658+55 
and BF1821+59, yield significant detections of power 
on short baselines. The flat-band power in these 
fields, determined from a joint analysis of the 0.63 — 
1.2 kA data, is found to be Qfiat = 18.9+ 7.2 /xK at 
68% confidence. To determine if this signal is due 
to point source contamination, we have reanalyzed 
the data assuming Q — Qa{t/la) rather than flat- 
band power. The effective multipole to which 
Qo is referred, is chosen so that Qo = Qfiat when 
both analyses are applied to the 0.63 — 1.2 kA data. 
We have determined Qo f° r several fields with known 
point sources and find that, for the long baselines 



(1.2 — 6.0 kA), the value of Qo is identical to that 
found from the short baselines (0.63 — 1.2 kA) used 
in the anisotropy analysis. Applying this analysis 
to the combined long baseline data for the fields 
BF0658+55 and BF1821+59, we find no evidence for 
point source emission and constrain Qo < 8.8 /iK and 
Qo < 15.8 fiK at 68% and 95% confidence. There- 
fore, we conclude with 95% confidence that Pois- 
son distributed point sources cannot be responsible 
for excess power found in the analysis of the short 
baseline data. However, if the signal is produced by 
several weak clustered point sources, such as in the 
PC1643+46 field, the scaling Qfiat oc ^ is not ob- 
served in either the real or simulated data. There- 
fore, weak clustered sources could be responsible for 
the observed signal. 
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Fig. 3. — The relative likelihood that the observed 
signal in each of the 1998 data sets is described by flat 
band power with amplitude Qfiat- The solid, dashed, 
dotted, and dot-dashed lines are for the BF 0028+28, 
HDF 1236+62, BF 1821+59, and BF 0658+55 data 
respectively. 

6. Expected Signals 

To interpret our results, it is informative to con- 
sider the level of the expected signals. In this section, 
we give our best estimates of the expected contribu- 
tions from primary and secondary CMB anisotropics 
and foreground sources to the observed Qfiat- These 
results are compiled in Table fj]. 

6.1. Primary Anisotropies 

We have convolved model CMB anisotropy power 
spectra with the window function of the experiment 
to determine the inferred flat band power signal, 



Q}lat 



\24tt J 



(16) 



The CMB power spectra were generated usi ng the 
CMBFAST code flSeljak fc Zaldarriaga 199^ ). The 



expected signal for this experiment is largely deter- 
mined by the total energy density of the universe, Oo- 
We have calculated Q fi a t for a range of values for fio 
while keeping the baryon density Qb — 0.05 and the 
Hubble constant h = 0.65 constant. For these ob- 
servations, the majority of the signal comes from the 
small region of overlap between the BIMA array di- 
agonal window function and the damping tail of the 
CMB. For total energy density £7q = 1-0 and 0.3, we 
expect primary anisotropy signals of Qfiat = 1.1 and 
5.8 respectively. 



Fig. 4. — The relative joint likelihood that the data 
from each year is described by flat band power with 
amplitude Qfiat- The dotted, dashed, and solid lines 
correspond to the 1997, 1998, and combination of 
1997 & 1998 data respectively. 



6.2. Secondary Anisotropies 

Sometime after recombination at redshift z ~ 
1100, the universe was reionized. We know that this 
ionization was essentially complete by redshift z ~ 5 
because spectra of distant quasars do not show a con- 
tinuum of absorption by neutral hydrogen (Gunn & 
Peterson 1965). The interaction of the CMB with 



the reionized universe leads to secondary anisotropies. 
There are three types of secondary anisotropies which 
are expected to make significant contributions on ar- 
cminute scales: the Vishniac effect; inhomogeneous 
reionization; and the SZ effect. 

The Vishniac effect is a "second order" Doppler 
shift which produces CMB temperature anisotropy 
by the fact that the large scale velocity field is mod- 
ulated by small-scale variations in baryon density 



( |Vishniac 1987| ; pstrikcr fc Vishniac 1986[ ). While the 
other secondary anisotropies discussed in this section 
require variations in ionization, the Vishniac effect 
acts in a universe that is uniformly ionized. Hu & 
White (1996) have determined the size of the effect 
for a range of reionization histories in the context 
of a critical CDM model. For reionization epochs 
z r = 5 and 100, they find signals of amplitude 
Qfiat ~ 1.7 and 3.6 /iK which peak at multipole mo- 
ments I ~ 5000 and 10000, respectively. 

Inhomogeneous reionization will imprint Doppler 
shifts, due to the velocities of the reioniz ed regions, o n 
the Compton scattered CMB photons (Kaiser 1984). 
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3000 4000 6000 8000 10" 
Multipole 1 

Fig. 5. — The "data weighted" window functions Wg. 
The dotted, dashed and solid lines correspond to the 
1997, 1998, and combination of 1997 & 1998 data 
respectively. 

The effect of inhomogeneous reionization has been 
studied most recently by Gruzinov & Hu (1998) and 
Knox, Scoccimarro, & Dodelson (1998). There are 
considerable uncertainties in the details of the gen- 
eration of ionized regions. To obtain an accurate 
result, the correlation of the ionizing regions must 
be taken into account. For universes that reionizes 
at zi — 26 and 31 , Knox, Scoccimarro, & Dodel- 
son (1998) find a flat-band temperature anisotropy of 
Qfiat ~ 1-8 and 2.5 /zK at £ = 5500. Observation of 
this signal would provide useful constraints on what 
are presently highly speculative reionization scenar- 
ios. 

The majority of luminous matter in massive clus- 
ters of galaxies is observed to exist in the form of 
ionized gas which has been heated by gravitational 
infall. This hot gas can present a considerable in- 
verse Compton scattering cross section to CMB pho- 
tons. The resulting spectral distortion in the direc- 
tion of a cluster of gala xies is known as the Sunyaev - 
Zel'dovich effect (SZE) (|Sunyaev fc Zel'dovich 19721) . 



The change in RJ CMB temperature in the direc- 
tion of a massive cluster can be as large as 1 mK. 
For a recent review see Birkinshaw (1999) and the 
references within. Several authors have computed 
the expected CMB anisotropy power spectrum due 
to the SZE in clusters of galaxies ( Atrio-Barandela fe 
Miil cket 199"9] [Komatsu fc Kitayama 1999| ; |Holder fc 
Ca dstrom 1999] ). The treatments differ in the range 
of cosmological models considered and the models for 
the cluster evolution. One general conclusion is that, 



at the small angular scales relevant for this experi- 
ment, the majority of the signal is due to distant less 
massive clusters and removing either the bright SZE 
or x-ray sources does not appreciably change the re- 
sults. In general, the results depend sensitively on 
the assumed cosmology and cluster evolution model. 
For example, Holder & Carlstrom (1999) find that at 
I ~ 5000, Q/iat — 1.3 — 8.0 ^K for the range of models 
they consider. 

6.3. Undetected Radio Point Sources 



In Section BTl we described our mechanism for 
measuring and removing point sources. As seen in 
Tables || and |[ the flux sensitivity of the data with 
baselines longer than 2.4 kA is comparable to that of 
the 0.63 — 1.2 kA data. Although we see no evidence 
for additional point sources, we cannot reliably con- 
strain the presence of point sources below 3cr in the 
high resolution maps which corresponds to a flux of 
- 300 - 500 fiJy. 

We have attempted to quantify the expected sig- 
nals from point sources both analytically and through 
simulations. The integrated source counts with fluxes 
less than S cu t have been measured at 8.4 GHz by Par- 
tridge et al. (1997). We take their result and scale 
it to the BIMA observing frequency of 28.5 GHz by 
using the average measured radio power law index 
(a — 0.77) from Cooray et al. (1998). The number 
density of sources then becomes 



N(> Sc„ t ) 



20 



aremm (s. 4GH 



-0 



(17) 

where (3 = 1.2. We follow the treatment of Scott & 
White (1999) and estimate the contribution to the 
power spectrum to be 



N(> S cut )S 2 c 



(18) 



where S cu t is the minimum source flux which we can 
remove from our maps. Using equation [l(], we can 
then determine the contribution of point sources to 
Q fiat ■ For the maximum residual source flux, S c ut — 
400 /iJy, we find Q/iat = 6.6 /iK. 

In order to test this approximation, we simulated 
distributions of point sources on the sky. For each 
simulated sky, we generated a sample of point sources 
with fluxes less than S cu t = 400 fi3y drawn from the 
dN/dS distribution given in Partridge et al. (1997). 
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The sources are placed at random in the field. These 
model skies are then Fourier transformed and added 
to a unique manifestation of visibilities consistent 
with the weights of one of our complete data sets. The 
simulations are then analyzed exactly as the real data. 
When we analyzed 100 manifestations of the sky gen- 
erated in this way, we found Qfiat = at 
68% confidence. So, the source simulations and ana- 
lytic approximation predict similar signals which, in- 
terestingly, are of the same order as the most likely 
signal in the data. 

However, as described in Section [l], the observed 
fields were selected to be free of bright radio sources 
at 1.4 GHz. If the sources have a falling spectrum, 
we will have selected fields with significantly lower 
than typical point source confusion at 28.5 GHz. Also, 
for the fields PC1643+46 and VLA1312+32, addi- 
tional OVRO observations were used to remove point 
sources down to ~ 200/iJy. Therefore, the estimates 
for point source contamination presented here are up- 
per limits to the expected signal in our data. 

6.4. Anomalous Dust Emission 



Recently, anomalous foreground emission at mi- 
crowave frequencies has been observed which is found 
to be strongly correlated with IRAS 100 /im m aps 
( |Lcitch ct al. 19971 ; |de Oliveira-Costa et al. 1997| ). It 
has been proposed that this emission may be due to 
either free- free emission ( Kogut et al. 1996 ) or dipolc 
emission from rapidly spinning dust grains ( Draine & 
Lazarian 1997). From a compilation of experimen- 
tal results, Kogut (1999) has determined a scaling 
between emission in the IRAS 100 /im band and at 
microwave frequencies. At 28.5 GHz, we expect this 
scaling to be approximately 17/xK/(MJysr _1 ). As 
mentioned in Section 2.2, we selected fields to have 



minimal 100 /im emission and contrast. We have de- 
termined the rms 100 /im flux for each of our observed 
fields. The resolution of the IRAS maps is 1.5' and 
therefore well matched to the angular scale on which 
the BIMA experiment is sensitive. The observed 
fields are found to have a range of rms fluctuations 
A/iotv = 0.04 - O.OgMJysr" 1 . Therefore, we ex- 
pect a rms temperature signal from this foreground of 
AT < 1.7/xK, which corresponds to Qfiat < 1-1 A*K- 

6.5. Systematic Errors 

The detected signals could also be the result of sub- 
tle systematic errors. The observations presented here 



represent the deepest images we have made of fields 
without strong SZE decrements due to known galaxy 
clusters and therefore could be subject to undiscov- 
ered systematic errors. Without success, we have ex- 
tensively searched for a non-astronomical explanation 
of the observed excess power. The results are found to 
be constant across the observing frequency band, in- 
dependent of baseline or telescope, reproducible from 
day to day, and uncorrelated with the position of the 
sun or moon during our observations. If this work is 
subject to systematic errors, deeper observations will 
be necessary in order for them to manifest themselves 
in a significant manner. 

Expected Contributions to Q fi at 



Signal 



Primary Anisotropy 


1.1 - 


5.8 


Vishniac Effect 


1.7- 


3.6 


Inhomogeneous Reionization 


1.8- 


2.5 


Sunyaev-Zel'dovich Effect 


1.3- 


8.0 


Point Sources 


< 6.6 


Spinning Dust/Free-Free 


< 1.1 


Total 


3.0- 


12.7 



Table 7: Expected flat-band power due to primary 
anisotropy, secondary anisotropy, and foreground 
confusion. 

7. Comparison with Previous Work 

There here have been many previous searches for 
anisotropy in the CMB at arcminute scales. The re- 
sults of this earlier work have been expressed in sev- 
eral different ways. Until recently, it was common for 
experimenters to quote limits on CMB anisotropies 
with a Gaussian autocorrelation function (GACF) 



C(9) = C exp{- — 



(19) 



where 6 C is the coherence angle and ^/Cq is the vari- 
ance of the CMB. Given the diagonal elements of 
the average window function, it is simple to convert 
betw een flat band power and GACF results ( Bond| 
1995| ). Here we express our results in terms of limits 
on temperature anisotropy with a GACF in order to 
facilitate comparison with the results of other experi- 
ments. At the scale of maximum sensitivity 9 C — 0.9', 
our 68% and 95% confidence limits on \/Ca/T cm b are 
6.5 x 10~ 6 and 9.6 x 10~ 6 respectively. 
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In previous work, both single dish and interfero- 
metric techniques have been used to perform sensitive 
searches for CMB anisotropics on arcminute scales. 



In Table B, we list the frequency, sky coverage, co- 



Compact Array (ATCA) has recently produced what 
were, previous to this work, the lowest limits on ar- 
cminute scale CMB anisotropies (Subrahmanyan et 
al. 1998). They observed at a lower frequency of 



hcrcnce angle corresponding to maximum sensitivity, 
and 95% confidence limits on variance and flat-band 
power for each of the most sensitive experiments and 
compare them with our results. 

Operating at a frequency of 20 GHz, the Owens 
Valley Radio Observatory (OVRO) 40 m dish has 
been used to measure sensitive differences between 
beams of ~ 1.8' fwhm separated by ~ 7' (Readhead et 
al. L989). They express their results in terms of limits 



on fluctuations with a GACF. At the coherence an- 
gle for which the experiment is maximally sensitive, 
9 C = 2.6', they constrain \/C^/T cmb < 1.7 x 10~ 5 
at 95% confidence. The effective total sky coverage 
of the experiment is estimated to be ~ 60arcmin 2 . 
More recently, the OVRO Ring experiment used the 
OVRO 40 m telescope to make a significant detection 
of anisotropy in a fi eld near the North Celestial Pole 
( Myers et al. 1993 ). These results are inconsistent 
with the earlier work at OVRO and are likely to be 
the result of foreground contamination. 

The SuZIE experiment is a drift scanning bolome- 
ter array. Fielded at the CalTech Submillimeter Ob- 
servatory (CSO), it was used to map ~ 213arcmin 2 
of blank sky at its operating frequency of 142 GHz 
(Church et al. 1997). Unlike the other experiments 
discussed here, radio point sources are not a signifi- 
cant source of confusion for SuZIE. They also express 
their results in terms of limits on fluctuations with a 
GACF. At the coherence angle of maximum sensitiv- 
ity, 6 C = 1.1', they find VO)/T cm6 < 2.1 x 10~ 5 at 
95% confidence. 

Interferometers have also proved to be very ef- 
fective for making sensitive maps of the sky with 
arcminute resolution. Using the Very Large Array 
(VLA) at 8.4 GHz, Partridge et al. (1997) obtained an 
extremely deep 21 arcmin 2 image of the sky. On the 
scale at which the experiment is most sensitive (reso- 
lution ~ 60") they find AT/T cmb < 2.0 x 10~ 5 at 95% 
confidence which corresponds to Qfiat < 35.2 [iK. In 
order to achieve this limit, they are forced to sub- 
tract a statistical estimate of the image variance due 
to point sources. Because we do not have exact knowl- 
edge of the window function of the VLA observations, 
we cannot determine the response of their system to 
CMB anisotropy with a GACF. 

A group working with the Australian Telescope 



8.45 GHz, but the larger size of the ATCA dishes com- 
pensates to make the window function of the ATCA 
and BIMA experiments similar. Using a single deep 
pointing of their array (~ 28 arcmin 2 ), they con- 
strained Qfiat < 23.6 [iK at 95% confidence on an 
angular scale corresponding to £ e // ~ 4600. They 
also express their results in terms of anisotropy with 
GACF. At the coherence angle for which the exper- 
iment is maximally sensitive, C = 1.0', they find 
VCo/Tcmb < 1-6 x 10" 5 at 95% confidence. 

8. Conclusion 

We have used the BIMA array in a compact con- 
figuration at 28.5 GHz to search for CMB anisotropy 
in seven independent fields. With these observa- 
tions, we have placed the lowest limits on arcminute 
scale CMB anisotropies to date. These results are 
determined from ~ 240 arcmin 2 of sky; this is the 
largest sky coverage of any of the arcminute scale 
anisotropy experiments. In the context of an assumed 
flat band power model for the CMB power spectrum, 
we find Qfiat = 5.6+ 5 ;g//K at 68.3% confidence and 
Qfiat < 14.1 /zK at 95.4% confidence with sensitiv- 
ity centered about harmonic multipole £ e f / = 5470. 
This result includes the three fields observed in 1997 
which were previously suspected to contain possible 
distant clusters. None of these fields contribute to the 
observed excess power. 

A detection of excess power is not surprising when 
one considers the level of the signals expected from 
secondary CMB anisotropies and foreground emis- 
sion. We have ruled out Poisson distributed point 
sources as the cause of the detected excess power 
at greater than 95% confidence, although, faint clus- 
tered sources could still be responsible. It is possible 
that we have detected some combination of secondary 
CMB anisotropies and faint radio point sources, how- 
ever, the confidence of the detection is only 44%. 

In the coming year, we plan to expand our ob- 
servations to include greater sky coverage and deep 
searches for point sources. Future observations with 
broader correlated bandwidth could reach sensitiv- 
ities an order of magnitude higher than presented 
here. With proper characterization of foregrounds, 
these observations may be able to place interesting 
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Comparison With Previous Work 



Experiment v (GHz) £l s ky (arcmin 2 ) 6 C (arcmin) 



95% Confidence Limits 

VCo/Tcmb Qflat 



2.1 x ict 5 

1.7 x 1(T 5 

1.6 x 10~ 5 
9.6 x 1(T 6 



SuZIE 
OVRO 40 m 
VLA 
ATCA 
BIMA 



142 
20 
8.4 
8.4 
28.5 



213 
- 60 
20 
28 

240 



1.1 
2.6 
1.0 
1.0 
0.9 



35.2 
23.6 
14.1 



Table 8: Frequency, sky coverage, coherence angle, and 95% confidence limits on the variance and flat-band power 
from previous work and the BIMA results. 



constraints on models for the reionization of the uni- 
verse. 
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